Observation of spin-wave characteristics in the two-dimensional ferromagnetic 

ordering of in-plane spins 
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The role of dipolar interactions and anisotropy are important to obtain, otherwise forbidden, ferro- 
magnetic ordering at finite temperature for ions arranged in two-dimensional (2D) arrays (monolay- 
ers). Here we demonstrate that conventional low temperature magnetometry and polarized neutron 
scattering measurements can be performed to study ferromagnetic ordering of in-plane spins in 2D 
systems using a multilayer stack of non-interacting monolayers of gadolinium ions. The spontaneous 
magnetization is absent in the heterogenous magnetic phase observed here and the saturation value 
of the net magnetization was found to depend on the applied magnetic field. The net magnetiza- 
tion rises exponentially with lowering temperature and then reaches saturation following a Tln(/3T) 
dependence. These findings verify predictions of the spin-wave theory of 2D in-plane spin system 
with ferromagnetic interaction and will initiate further theoretical development. 

PACS numbers: 75.70.Ak, 75.60.Ej, 75.50.Xx 
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Ferromagnetic materials confined in ultra-thin films 
and multilayered structures are being studied extensively 
for the development of high-density magnetic data stor- 
age devices and to refine our basic knowledge in low- 
dimensional physics P, 0, IE lit- Recent advances in 
growth techniques such as molecular beam epitaxy and 
magnetization (M) measurement techniques based on the 
magneto-optical Kerr effect have enabled us to measure 
small magnetic signals as a function of applied magnetic 
field (H) and temperature (T) even from one atomic 
monolayer of a ferromagnetic material and a wide range 
of ordering effects has been observed 0, IE IE IS EH • 
These measurements have also demonstrated the exis- 
tence of a spontaneous magnetization and have revealed 
hysteresis curves in two- U and one-dimensional sys- 
tems, where magnetic ions are arranged in a grid or 
in a line within a monolayer. A recently generalized 
theorem [E showed, following spin- wave theory, that 
long-range ferromagnetic order and hence spontaneous 
magnetization cannot exist at finite temperature in a 
two-dimensional systems provided spin-spin interactions 
are isotropic and short range. A theoretical formalism 
EE IE EH 113 and computer simulations [E have been 
developed to include anisotropy and dipolar interactions 
to explain the apparent contradiction between theory and 
experiment in low-dimensional systems. A 2D array of 
magnetic ions with lattice parameter 'a' of spins S can 
be described by a Hamiltonian, 



these terms in equivalent magnetic field units as, 



= JS, 2/iFjHd = iiragS, 2/j, B H k = 6KS (2) 

In the above expression a(~ 1) depends on the lat- 
tice type and g is equal p to (2/ie) 2 /a 3 . K is the 
anisotropy constant. The magnetization reduction due 
to thermally activated spin waves was calculated with 
this Hamiltonian and one obtains a non-zero temperature 
for long-range ferromagnetic ordering as a gap of width 
A 2 = 2/i_B_ff£ ff opens up at the bottom of the spin wave 
spectrum for an easy-magnetization axis (z) perpendic- 
ular to the film plane. The easy-magnetization axis is 
determined by the sign of the effective anisotropy field 
(Hf = H k - H d ) defined Q by 
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(6K oS S) with K eS = K- 



lixtxg 



(3) 



H = Hex + ~Hd + T~Lk 



(1) 



This explains the observation of hysteresis curves in a 
monolayer with spins oriented normal to the surface 
[E IE m- However, the situation is quite different for 
spins oriented in an in-plane direction with < 

as the spin- wave spectrum remains gapless. The long- 
range character of the dipole interactions was found 
to be responsible for creating a pseudo-gap A xy = 
(irSg/2)^y (6\K c s |/J) in the spin-wave spectrum that 
may give rise to long-range ferromagnetic order in 2D in- 
plane spins provided \K c s\ > K c = n 2 g 2 /(6J). The tem- 
perature dependence of the magnetization M(T) above 
a transition temperature T c (— 6S\K c s\/ Kb) takes the 
same form 



The strength of the three terms arise from exchange, 
dipolar and magneto-crystalline anisotropic interactions 
respectively, and have been approximated by expressing 



M(T) = M [l - AT1n(0T)] 



(4) 



for the ordering of in-plane spins (with (3 = Kb/ £± X y = 
2y/Q\K c s | J I (jrgT c )) and out-of-plane spins (with (3 — 



2 



K B /A Z = l/T c ). Here A = K B /(4nJS 2 ) and M is 
the saturation value of the net magnetization that de- 
pends on the field applied to carry out measurements 
Below T c spin wave theory predicts |7j an enhance- 
ment of M(T) as [1 - A exp(-l/(/3T))] and [1 - CT y ] 
for out-of-plane and in-plane orderingrespectively where 
C depends on A xy and v expected ja, 0, EI to be 3/2. 
For < |-Keff| < 7I ' 2 5 2 /(6J) in-plane spins can not stabi- 
lize in a homogeneous phase as the magneto-crystalline 
anisotropy becomes large enough to pull some of the spins 
in the out-of-plane direction and create a ripple like in- 
stability 0,0- It is known that both the average magne- 
tization M(T) as well as the initial susceptibility x(T) is 
proportional to the physical extent of the ordered phase 
I* that minimizes the zero-field energy 6] and can be 
written as 

M oc HI* and \ « V with I* oc exp(~7T) (5) 

It is expected that at low enough temperature, I* reaches 
saturation either because I* becomes comparable to the 
sample size or due to a freezing of the ripple walls. The 
net magnetization M(T) of the ordered domains then 
should follow the spin- wave prediction (Eq. (4)) to reach 
saturation. 

Here we present the verification of this theoreti- 
cal prediction of 2D ferromagnetic ordering of in-plane 
spins using Langmuir-Blodgett (LB) films of gadolinium 
stearate. The presence of a large multilayer stack of non- 
interacting monolayers of gadolinium has enabled us to 
carry out conventional quantitative magnetization mea- 
surements at sub-Kelvin temperatures. We could also 
use polarized neutron scattering measurements [l^L fl"3 | , 
to show that the ordered phase of the in-plane spins is 
inhomogeneous and that the monolayers remain uncor- 
rected even when the net magnetization reaches a satu- 
ration value. 

In the metal-organic structure formed by LB tech- 
niques 0, 0, 0] , gadolinium ions are separated by ap- 
proximately 5A within a monolayer to form a distorted 
hexagonal 2D lattice and the monolayers are separated 
from each other by 49A by organic chains (Fig. 1). Films 
having 50 such monolayers of gadolinium ions were de- 
posited on 1 mm thick Si(001) substrates and character- 
ized using x-ray reflectivity technique, as discussed earlier 

Neutron reflection measurements (refer Fig. 2) were 
carried out on the CRISP reflectometer at the Rutherford 
Appleton Laboratory (RAD, UK using a cold, polychro- 
matic neutron beam [lj, and at the ADAM beamline 
[20I |2]| of the Institut Laue-Langevin (ILL), Grenoble, 
France, using a monochromatic cold neutron beam. In 
Fig. 2(b) we have shown spin-polarized neutron reflectiv- 
ity data taken at 2K by applying a magnetic field of 13 
kOe along an in-plane direction (refer Fig. 1). It is known 
that, if the polarization of the neutrons defined to be par- 
allel (+) or antiparallel (-) to the applied field direction 



(along +y axis) dEHH, the intensity of a multilayer 
Bragg peak in this geometry increases with effective scat- 
tering length b cff = b ± A/j, y with A = 0.2695 x 10~ 4 
A//is. In the left inset of Fig. 2(b) we have shown inten- 
sity profiles of parallel (+) and anti parallel (-) incident 
neutrons at the first peak position. The average of (+) 
and (-) profiles represent the non-magnetic contribution 
to the reflectivity. Systematic analysis of all these pro- 
files provide us the value of [i y , the component of aver- 
age moment per gadolinium ion along +y direction. The 
obtained values of fj>y CIS Ct function of H at 4.2K and 
at 1.75K using the CRISP and ADAM spectrometers 
respectively are shown in Fig. 2(a) along with data ob- 
tained from the magnetization measurements [ljj at 2K 
and 5K. Results of these two independent measurements 
show that the obtained average saturated moment per 
gadolinium ion is much less than the expected 7^b and 
confirms the existence of a heterogeneous phase. In the 
right inset of Fig. 2(b) we have shown the transverse dif- 
fuse neutron scattering intensity profile at the first Bragg 
peak. The hyper-geometric line shape profile confirms 
that the in-plane correlation is logarithmic in nature and 
that the interfaces are conformal |U 0| • It should be 
noted that unlike in x-ray measurements the scattering 
here originates primarily from the metal heads. The line 
shape and the associated parameters were found to be 
independent of T, H and hence exhibit the absence of 
conformality in the magnetic correlations between inter- 
faces. This again confirms that the LB films represent a 
collection of isolated 2D spin-membranes of gadolinium 
ions . 

Now we present the results of conventional magneti- 
zation measurements carried out at sub-Kelvin temper- 
ature to understand the nature of the ordering. These 
measurements were performed in a conventional way by 
measuring forces exerted on a sample situated in a spa- 
tially varying magnetic field with a Faraday balance [25j . 
In Fig. 2(a) we have shown M vs. H data taken at lOOmK 
and 500mK temperature in an in-plane (+y) direction. 
This data reconfirms the absence of hysteresis and re- 
manence (M = at H = 0). The saturation value of 
the net magnetization at lOOmK and 500mK found to be 
12.7 x 10~ 6 emu/mm 2 = 5.4/i^/ Gd atom; much lower 
than the expected 7.0 ^/Gd atom for a homogeneous 
ferromagnetic phase. 

In Fig. 3(a) we have shown the magnetization data 
taken with different fields as a function of tempera- 
ture. At higher temperatures these data were fitted 
with Eq. (5), and the expected exponential dependence is 
also observed in the magnetization data extracted from 
our neutron reflectivity measurements (inset of Fig. 3(a)). 
The values of 7 obtained from fitting were found to in- 
crease with the reduction of H and at 0.25 kOe it is found 
to be 2.162 K . It is also observed that at a lower field, 
the magnetization at a fixed temperature is nearly pro- 
portional to the applied field (5.03 x 10~ 7 at 0.25 kOe 
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and 1.29 x 10 -6 emu/mm 2 at 0.5 kOe and at a tem- 
perature 0.9 K) as predicted for inhomogeneous striped 
phases (refer Eq. (5)). The amount of majority phase 
grows exponentially as we lower the temperature for each 
applied field but below a certain temperature this growth 
stops as the walls of the majority phase freeze. Below 
this temperature measured data do not follow Eq. (5) 
and the net magnetization of the majority phase M(T) 
increases with lowering temperature following the ther- 
mally activated spin waves as given in Eq. (4). We fit- 
ted all the data with Eq. (4) and obtained Mq values 
as 0.9 x 10~ 6 , 1.6 x 10~ 6 ,3.2 x 10~ 6 ,4.5 x 10~ 6 ,7.9 x 
10~ 6 ,9.5 x 10~ 6 emu/mm 2 with 0.15, 0.25, 0.5, 1.0, 2.5, 
5.0 kOe magnetic fields respectively. These saturation 
values of net magnetization indicate that the percent- 
age of the ferromagnetic phase is increasing from 7.1% 
to 74.8% as we approach the maximum saturation value 
of net magnetization obtained of 12.7 x 10~ 6 emu/mni2 
(= 5.4 fiB /Gd atom) as shown in Fig. 2(a). We extracted 
the value of exchange J as 8.76 x 10 -19 erg (or H ex = 
0.165 kOe) from the fitted value of A{= 1.02 K' 1 ) for 
the 0.25 kOe data. We obtained (3 as 3.4 for all the data 
and hence |-Keff| was calculated to be 1.7 x 10~ 19 erg 
(or Hf = 0.19 kOe) for 0.25 kOe data giving T c =26 
mK. In this calculation g was 6.88 x 10 -18 erg, assum- 
ing that one gadolinium atom occupies 2.5 Ax20 A 2 , as 
obtained from neutron and x-ray analysis (refer Fig. 1). 
This proves that we are dealing with an inhomogeneous 
phase as < \K eS \ < K c {= 8.89 x 10~ 17 erg). Al- 
though Eq. (4) describes the temperature dependence of 
the magnetization for ferromagnetic ordering of both in- 
plane and out-of-plane spins, the argument of the loga- 
rithmic function can become less than 1 only for in-plane 
ordering. Unusually low values of H ex and H'jf with a 
rather large value of Hd(= 16.3 kOe) make [3T < 1 even 
for T > T c - this situation has not been reported ear- 
lier to the best of our knowledge. It is interesting to 
note that all the magnetization data shown in Fig. 3(a) 
attains respected saturation values Mq at temperature 
To = l//3(~ 0.29 K) and a maximum magnetization at 
temperature T m — l/(e/3) (~ 0.108 K). Our experimen- 
tal uncertainties below 100 mK prohibit us from com- 
menting on any reduction of magnetization below this T m 
but Eq. (4) with (3T < 1, all the data is fitted quite well. 
Further theoretical development is required to under- 
stand the thermal activation of spin- waves with in-plane 
ordering especially as we approach below T c (= 26 mK). 

In Fig. 3(b) we have shown zero-field-cooled (ZFC) and 
field-cooled (FC) magnetization data taken with 0.15 kOe 
and 0.5 kOe field along with the fitted curve from Eq. (4). 
We observe a blocking temperature (T&) of 125 mK be- 
low which there is branching in both the ZFC and FC 
data. This result reconfirms that the observed ferromag- 
netic ordering requires an external field to stabilize and 
such a low T]y indicates the existence of very small spin 
clusters in the ZFC phase. Application of the field lowers 



the activation energy required for the randomly oriented 
domains to increase the number of spins in the ferro- 
magnetically ordered majority phase and increases the 
activation energy of the reverse transition 6}. As a re- 
sult we do not observe Tj, in the measured temperature 
range for 0.5 kOe. 

In conclusion, we have demonstrated that polarized 
neutron scattering and conventional magnetization mea- 
surements can be used to study 2D ferromagnetic order- 
ing of in-plane spins using a stack of magnetically uncor- 
rected spin-membranes formed with gadolinium stearate 
LB film. The in-plane ordering observed here shows that 
a spontaneous magnetization is absent even at 100 mK 
and saturation value of the net magnetization increases 
with a lowering in temperature. The magnetization is 
found to increase exponentially with a lowering in tem- 
perature due to the exponential increase of the physical 
extent of the ferromagnetic domains in the heterogeneous 
phase. The ferromagnetic domains ultimately saturate 
following T \n((3T): characteristic of thermally activated 
spin- waves and are found to be valid for even [3T < 1. 
We believe these experimental results will initiate further 
theoretical development. 
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Figure captions: 

FIG. 1: (a) Schematic diagram of the out-of-plane and in- 
plane structure of the gadolinium stearate Langmuir Blodgett 
film is shown with the scattering geometry employed for the 
polarized neutron reflectivity measurements, x — z plane is 
the scattering plane and the magnetic field is applied along 
the +y direction. A and a are the wavelength of the radiation 
and angle of incidence respectively. 

FIG. 2: (a) In-plane Magnetization curves obtained as a func- 
tion of the field (H) using neutron reflectivity measured at 
4.2K (diamond) and 1.75K (star) compared with conventional 
magnetization data ll7|l measured at 2K (down-triangle) and 
5K (up-triangle). Solid lines are the fits with a modified Bril- 
louin function |l7l| . Magnetization measured at lOOmK and 
500mK is shown for the 1st (symbols) and 2nd (line) cycle 
of the hysteresis loop, (b) Neutron reflectivity data (sym- 
bols) at H = 13 kOe and at T—2K for the neutron spin 
along (+) and opposite (-) to the magnetic field direction 
with the corresponding fit (line). In the left inset the first 
Bragg peak is shown in (+) and (-) channels in an expanded 
scale. Right inset: Transverse diffuse neutron scattering pro- 
files (symbols) measured at 2K with unpolarized and polar- 
ized neutron beams. The solid line is a fitted hypergeometric 
curve as described in the text. 

FIG. 3: (a) Sub-Kelvin magnetization results with various 
applied fields (symbols) fitted with Eq. (4) (Black line) and 
with Eq. (5) (wine coloured dashed lines). Dotted lines indi- 
cate the temperatures T m and To (refer to text). Inset shows 
the magnetization obtained from neutron measurements as a 
function of temperature (symbols) and fit with Eq. (5) (line) . 
(b) ZFC (green circles) and FC (blue stars) along with the fit 
(line) for FC measurements. 
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